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Investigating the properties of AGN feedback in hot 
atmospheres triggered by cooling-induced gravitational 
collapse 
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ABSTRACT 

Radiative cooling may plausibly cause hot gas in the centre of a massive galaxy, or 
galaxy cluster, to become gravitationally unstable. The subsequent collapse of this 
gas on a dynamical timescale can provide an abundant source of fuel for AGN heating 
and star formation. Thus, this mechanism provides a way to link the AGN accretion 
rate to the global properties of an ambient cooling flow, but without the implicit as- 
sumption that the accreted material must have flowed onto the black hole from 10s of 
kiloparsecs away. It is shown that a fuelling mechanism of this sort naturally leads to 
a close balance between AGN heating and the radiative cooling rate of the hot, X-ray 
emitting halo. Furthermore, AGN powered by cooling-induced gravitational instability 
would exhibit characteristic duty cycles (5) which are redolent of recent observational 
findings: 8 oc Lx/cr^, where Lx is the X-ray luminosity of the hot atmosphere, and 
cr* is the central stellar velocity dispersion of the host galaxy. Combining this result 
with well-known scaling relations, we deduce a duty cycle for radio AGN in elliptical 
galaxies that is approximately oc Mbh 1 ' 5 , where A/bh is the central black hole mass. 
Outburst durations and Eddington ratios are also given. Based on the results of this 
study, we conclude that gravitational instability could provide an important mecha- 
nism for supplying fuel to AGN in massive galaxies and clusters, and warrants further 
investigation. 
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1 INTRODUCTION 

The cooling time of X-ray emitting gas near the centres of 
many massive galaxies and galaxy clusters is much shorter 
than the Hubble time. In the absence of heat sources, sig- 
nificant quantities of the gas would cool and form stars. 
However, X-ray spectroscopy has shown that the rate at 
which gas cools to low temperatures is significantly lower 



than first expected (e .g. Peterson et alj 2001 



2001; IXu et all [20021; Sakelliou etaiT2002: 



2003l ; lKaastra et al.ll2004l ; |Peterson fe Fabian 



Tamura et al 



Peterson et al 
20061 ) suggest- 



ing that the gas is somehow being reheated. 

Based on both observational and theoretical evidence, it 
is generally assumed that energy input by a central AGN is 
predominantly responsible for reheating the gas. For exam- 
ple, elliptical galaxies are commonly the hosts of powerful ra- 
dio Active Galactic Nuclei (AGN). These sources give rise to 
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lobes of radio emission embedded in the X-ray emitting gas 
whic h permeates massive galaxies and clusters of galaxies 
(e.g. iBirzan et al.ll20ol ; iBest et af1l2005l; IrJunn et al.1 [20051 : 
IRaffertv et ail 120061 : IBest et al.l 120071 : IShabala et all 120081 ). 
Moreover, recent observational studies of Brightest Clus- 
ter Galaxies (BCGs) suggest that radio AGN activity is 
related to the thermal state of its environment. Systems 
with short radiative cooling times, or a low central en- 
tropy, are more likely to exhibit active star for mation , opti- 
cal line-emission and j et-producing AGN (e.g . Burns!ll990l: 



Crawford et al. 19991: Cavagnolo et al.1 |20Q. 



Mitta l et all 



20091 ; IRaffertv et al.H2008l b This suggests that AGN activ- 



ity is part of a feedback loop that can prevent the ambient 
hot gas from cooling, and is likely to have other important 
consequences for its environment. 

Building o n the ideas of early w ork (e.g 



Binnev fc Taborl Il995l; [Tucker fc Daridl Il997l ; 
Ciotti & Ostrikcr 200l|; ISilk fc Reesl ll99ST ). theoretical 



studies have drawn attention to the potential, wide-ranging 
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impact of AGN feedback. For example, semi-analytic models 
of galaxy formation have demonstrated that, in principle, 
AGN heating can both reheat cooling flows and explain 
the exponential cutoff a t the bright end o f the galaxy 
luminosity function (e.g . iBenson et all 120031; ICroton et al.l 
120061 : iBower et all 12009 ). see also ijShort fc Thomas! l2009h . 



More recently, AGN heating has been shown to be crucial 
in shaping the X-ray luminosity-temperatu r e relation of 
mass i ve galaxies (e.g. iPuchwein et alj 120081 ; IBower et al.l 
l200Sl ; |Popell2009l ). 

Despite these findings, the fundamental details of AGN 
feedback remain poorly constrained. Most significantly, 
there is no clear consensus on how information about the 
thermal state of the X-ray emitting atmosphere is transmit- 
ted to the black hole at the centre of the galaxy, or clus- 
ter. Understanding this link is of great importance because 
it probably facilitates the observationally-inferred long-term 
balance between AGN heating and gas cooling, and almost 
certainly governs the duty cycle of AGN activity - that is, 
the fraction of time an AGN is active. 

Generally, AGN feedback is assumed to be powered by 
one of the following: a) Bondi accretion of hot material 



in the vicinity of the b l ack hole (e.g. Cattaneo fc Tev ssicr 
20071 ; ISiiackf et al. l2008l ; |Puchwein et a l. 2008; F abian et al 



20101); b) material directly from the ambient cooling flow 
(e.g.lPizzolato fc Sokeij|2005l ; IBower et al]|2008l ; |Popell2009l ; 



iPizzolato fc Sokerj l201oi ~ While there are plausible argu- 
ments for both, based on our current interpretation of the 
available observations, there are also difficulties. For exam- 
ple, the power output from Bondi accretion is unlikely to be 
sufficient to balance the gas cooling rates in massive clus- 
ters (e.g |Sokerll2006l ). The biggest problem with models in 
category b) is finding a mechanism by which the central 
black hole receives information from the cooling flow about 
the thermal state of the ICM. For example, it is difficult to 
conceive of how material might flow all the way from the 
cluster cooling radius onto a central black hole, such that it 
prompts an AGN outbu rst on a useful timescale (though see 
IPizzolato fc Soker|[2005l '). 

The model investigated here describes a mechanism 
that c an reconcile the problem f aced b y models of type b). 
While IPizzolato fc Sokerj |2005l . l201dh focused on thermal 
instability of the ICM as a mechanism for delivering fuel to 
the AGN, we focus on the characteristics of AGN feedback 
that is triggered when the hot gas which resides near the 
centre of a galaxy becomes gravitationally unstable. That 
is, gas which was previously self-supporting against gravity 
is somehow destabilised and falls freely towards the galaxy 
centre, on a dynami cal timescale , where it forms sta r s and 
fuels the AGN (e.g. ISilk fc Reesl Il998l ; iFabianl 120091 ; iKind 
2009). Grav itational instabiliti es can be induced by merger 
events (e.g. ISilk fc R ccs 199 8]) and t he effect of radiative 
cooling (e.g. iBirnboim fc Dekeill2003l ). Of these two possi- 
bilities, only instabilities induced by gas cooling can lead to 
self-regulated AGN heating, but it is important to mention 
that merger-driven AGN feedback will delay the onset of 
gravitational instability induced by radiative cooling. The 
extent to which this happens is unclear, but the effect may 
be more significant in group and field environments (e.g. 
iKavirai et alB oiO), than in clusters. 

Perhaps importantly, AGN heating driven by cooling- 
induced gravitational instabilities has the potential to be pe- 



riodic, as explained by the following argument. In massive 
galaxies and clusters there is a continual inflow of material 
from the large-scale environment, due to the effects of the 
ambient cooling flow. As the mass of material in the hot 
atmosphere near the centre of a galaxy builds up, it can be- 
come gravitationally unstable, meaning that some fraction 
of the gas will collapse and flow on a dynamical timescale 
to the centre of the gravitational potential, thereby fuelling 
the AGN. Once this fuel has been consumed, the AGN will 
be starved of new material until the nearby hot gas once 
again becomes gravitationally unstable. A second instabil- 
ity will commence after the external inflow has delivered a 
sufficient quantity of new material. Thus, for this mode of 
fuelling, the typical duration of an AGN outburst is largely 
governed by the local dynamical timescale, while the time 
between the onset of successive outbursts is controlled by the 
ambient cooling flow. This is notable because observations 
indicate that AGN are only active for a fraction of time de- 
term ined by their environment and host g alaxy properties 
(e.g. iBest et al.l 120051; IShabala et ai1l2008r i. In addition, as 
demonstrated in iPopel l|201li r periodic heating can supply 
the energy required to balance gas cooling, with the mini- 
mum effort. Thus, periodic AGN activity appears to be an 
energetically favourable heating strategy. Therefore, start- 
ing from the assumption that AGN fuelling is attributable 
to gravitational instabilities induced by gas cooling, we de- 
rive AGN duty cycles, outburst durations, star formation 
rates and discuss implications for numerical simulations. 

The outline of the article is as follows. In Section 2 
we present the main model and place constraints on the 
phenomenological parameters. In Section 3 we show the ex- 
pected duty cycles and outburst durations that would be as- 
sociated with AGN heating that is driven by cooling-induced 
gravitational instability. The implications are discussed in 
Section 4, and the results are summarised in Section 5. 



2 OUTLINE OF THE MODEL 

In this investigation, we study the properties of AGN feed- 
back that is driven by the gravitational collapse of hot gas 
in the vicinity of the supermassive black hole. Therefore, 
we must first describe the physical conditions of the gas 
in massive galaxies and clusters that can justifiably lead to 
gravitational instability. Following this we define the charac- 
teristic timescales of the system and the AGN power output 
expected from this mode of fuelling. These quantities are 
then used to determine the observable features predicted by 
this scenario. 



2.1 The onset of gravitational instability 

A self-gravitating sphere of gas will become gravitationally 
unstable when the sound-crossing time is greater than the 
gravitational free-fall time. This is the well-known Jeans cri- 
terion. Strictly speaking, the criterion only applies in the 
limit that the polytropic index of the gas - 7 e ff, as de- 
fined in equation (U — is also lower than som e critical value, 
7crit ~ 1-2, (e.g. see IBirnboim fc D ckcl 2003, and references 
therein). Since the actual adiabatic index of an ideal, non- 
relativistic monatomic gas is 7 = 5/3, a gas of this type 
can only become gravitationally unstable if its behaviour is 
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modified by heating, cooling and work. For example, if a 
parcel of gas is heated by an amount AQ, and then does an 
equal amount of work on its environment (i.e. AW = AQ), 
the gas is defined as isothermal, since it acts to maintain a 
constant temperature. In this case, the polytropic index of 
the gas must be y e g = 1, by definition. 

In the case of the hot, X-ray emitting gas that resides in 
massive galaxies and galaxy clusters, we envisage a central 
gas mass that is embedded in a pressurised ambient medium. 
Consequently, the scenario is somewhat reminiscent of the 
isothermal sphere which is prone to gravitational instability 
if its mass exceeds the Bonnor-Ebert limit. However, it is 
important to note that in this study we do not make any 
assumptions about the density distribution of the hot gas. 

Observations clearly indicate that this gas must be sub- 
ject to radiative cooling, heating from AGN and stars, and 
the work done by its surroundings. Therefore, it seems ex- 
tremely likely that 7 c ff =^5/3 for the X-ray emitting gas. 
However, it is unclear exactly what values the polytropic 
index might take because the processes acting on the gas 
are highly uncertain. Nevertheless, consider the behaviour 
of a gas mass near the centre of a massive galaxy clus- 
ter. As the gas radiates energy it loses pressure support; 
the weight of the overlaying gas then does work and com- 
presses it. The contraction can only persist if the compres- 
sional heating resulting from the inflow is radiated away on 
a timescale that is shorter than the dynamical timescale, i.e. 
icooi < idyn- This is perfectly possible because a gas parcel 
of temperature T and number density n, radiating via ther- 
mal bremsstrahlung, has a cooling time that is tx T 1 ^ 2 /n. 
As the gas cools, the temperature will fall, and the density 
will rise, so that the cooling time shortens. It is this con- 
dition (t coo i < tdyn) that explains observations which show 
that the temperature of the gas remains lower near the cen- 
tre of the cluster than further out. Indeed, it is precisely 
this argument that predicts a cooling catastrophe, and the 
need for AGN feedback. However, as shown below, the inflow 
condition also predicts the onset of a cooling-induced gravi- 
tational instability which provides a potentially informative 
description of AGN fuelling in hot atmospheres. 

Thermodynamically, the mass of gas described above 
behaves like a parcel that cools down (rather than heats 
up) when compressed. The appropriate polytropic relation 
linking the temperature, T, and volume, V, of the gas parcel 
is T\ ry " ft ~ 1 — constant. Thus, the hot gas in the centres of 
massive galaxies and clusters acts as if < j e ff < 1. Since 
this is less than all plausible values of the critical polytropic 
index, the hot gas is likely to be susceptible to gravitational 
instability Q. 

For completeness, we also present a brief derivation to il- 
lustrate the general phenomena that influence the polytropic 
index of a gas. Th i s der ivation closely follows the work of 
iBirnboim fc Dekel ([2003) who define the polytropic index, 
7 c ff , of a Lagrangian fluid element in terms of the logarith- 
mic time derivative of its pressure, P, and density, p, such 
that 



7eff 



dlnP/dt 
d In p/dt ' 



(1) 



1 The gas exhibits characteristics which are somewhere between 
isothermal (7 e g = 1) and isobaric (7 c ff = 0) states. 



Accounting for gas cooling, the polytropic index can be ex- 
pressed as 



Tefl = 7 



n 2 A(T) 
pe 



(2) 



where 7 is the actual adiabatic index of the gas, p is the rate 
of change of gas density and e is the internal energy per unit 
mass. As usual, n is the number density of the gas, T is the 
gas temperature and A(T) is the cooling function. 



2.2 Characteristic timescales and flow rates 

To calculate the local dynamical timescale, we refer back 
to the definition of the Jeans instability criterion. A self- 
gravitating gas mass will become gravitationally unstable 
when its mass exceeds some critical value that is compa- 
rable to the Jeans/Bonnor-Ebert masses. In the centre of 
a galaxy, the self-gravity of the gas will become important 
when its characteristic velocity dispersion is comparable to 
the central stellar velocity dispersion, cr*. Thus, if a gas mass 
M with radius R located at the centre of a galaxy becomes 
gravitationally unstable, it will collapse on the dynamical 
timescale, ta yn = R/<J*. The average mass flow rate during 
this time will be 



M 



AM 
~Af 



M 

^dyn 



R 



(3) 



We note that, for a self-gravitating cloud, the characteris- 
tic velocity dispersion of the constituent particles is related 
to the gravitational potential energy by the virial theorem: 
aGM/R — a 2 , where G is Newton's gravitational constant. 
In this description the numerical constant a ~ 1 depends on 
the density distribution of the gas. Thus, the mass flow rate 
can be expressed as 



M 



~e^G 



1000 



200 kms- 



M Q yr" 



(4) 



This is the well- known form of the dynamical mass flow rate 
(c.f. King 2009), modified slightly to account for arbitrary 
density distributions. From this it is straightforward to show 
that the duration of the collapse will be M/M — td yn - 

Following the collapse, mass will flow towards the centre 
of the galaxy's gravitational potential where a fraction is 
likely to be accreted by the supermassive black hole, and 
the remainder presumably forms stars nearby or is expelled 
from the galaxy as a result of feedback. 

As described in the introduction, the fuelling rate pro- 
vided by gravitational instability can be periodic, for the 
following reason. The postulated gas cloud near the galaxy 
centre acts as a reservoir; this reservoir is depleted when the 
cloud collapses, but is replenished by the inflow of new ma- 
terial from the larger-scale environment, which forms a new 
cloud. When it reaches the critical mass, the second cloud 
will also become gravitationally unstable and collapse, and 
so on. 

The provenance of the material that replenishes the col- 
lapsing gas cloud is not clear, except that it must come from 
the ambient hot atmosphere. The simplest possibility is that 
the cloud is built up during two phases: 1) an initial collapse 
of material from slightly further out in the hot atmosphere; 
2) the subsequent slow inflow of additional material due to 
the ambient cooling flow. If this is the case, the first phase 
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of growth of the new cloud will occur on a similar dynamical 
timescale to the gravitational collapse of the previous cloud. 
However, at these early times, the new cloud is likely to be 
hotter and of lower mass than the previous cloud, ensuring 
that it can be gravitationally stable. More precisely, as long 
as 7eft < 7crit ~ 1-2, and the cloud mass is below the crit- 
ical mass, the cloud will not collapse. The second phase of 
growth will occur on the timescale required for the cooling 
flow to build up the cloud mass to its gravitationally unsta- 
ble limit, at which point a collapse will be initiated. Using 
this argument, we can estimate the average time between 
gravitational collapses, as shown below. 

Overall, the local mass inflow rate, M ex t, from the cool- 
ing flow will be a slowly varying quantity governed by the 
difference between the time-averaged heating and cooling 
rates. Then, in the limit that the majority of the cloud mass 
is built up during this phase, the characteristic time between 
the triggering of successive gravitational collapses must tend 
to 

r=^. (5) 

Mext 

Assuming that the AGN is only fuelled while the cloud col- 
lapses, it will be active for a fraction of time S = tdyn/T, 
known as the duty cycle. By combining equations Q and 
(JS|, we obtain a simple functional form for the AGN duty 
cycle without having to explicitly calculate the AGN heating 
rate 



T 



CtGMe, 



Ma, 



(6) 



Thus, in the present model, the AGN duty cycle can be 
straightforwardly related to the local gravitational poten- 
tial, through a,, and the external environment through 
Mext - This is potentially significant, because oberv ations 
(e.g. iBest et al.l [20051 ; lBestll2007l ; IShabala et a"i1l2008l ) indi- 
cate that the radio AGN duty cycle is heavily influenced by 
both local and environmental effects. Below, we argue that 
Af ext is probably closely related to the mass inflow associ- 
ated with the ambient cooling flow. 



2.3 AGN heating rates 

If a fraction, f3, of the inflowing mass rate, M, reaches the 
black hole, the accretion power output will be 



H = rjpMc = rj——c 



10 4 ' ( ilSL ) I 11 ) ergs" 

V10 - V \200kms- V 6 

(7) 

where the assumed accretion efficiency is r\ ~ 0.1. The char- 
acteristic value of /3/a = 10 -3 is motivated below, and yields 
favourable comparisons with the observationally-inferred ra- 
dio AGN duty cycle. 

Clearly, the AGN energy injection rates associated with 
this fuelling mechanism can be very large. Such values are 
considerably larger than would be expected to arise from 
the accretion of nearby hot gas and, if observed in a real 
system, would more commonly be associated with merger- 
driven fuelling events. However, the reasoning above indi- 
cates that such an interpretation is not necessarily exclu- 
sive - high AGN fuelling rates can also be a consequence 
of gravitational instability resulting from gas cooling. In ad- 
dition, it is important to remember that values estimated 



from equation (0 represent the instantaneous heating rate 
- the time-averaged values are much more modest. 

By definition, the time-averaged AGN power output is 
written H = 8H, which can be expanded using equations 
© and © to give 

H = 8H = ril3M^ t c 2 . (8) 

Equation (|8]) shows that the time-averaged AGN heating 
rate depends on the properties of the large scale envi- 
ronment, through Afoxt- From this we also conclude that 
Mext will evolve until the time-averaged heating rate closely 
matches the ambient cooling rate of the gas. In the limit 
that the time-averaged AGN heating does balance gas cool- 
ing, H = Lx, equation ((8} shows that 



175 



10 45 ergs" 1 



10- 



M yr \ 



(9) 

where again we have assumed r\ w 0.1 and /3/a ~ 10 3 , with 
a ~ 1. Below, we motivate the constraints on j3 by consid- 
ering the Eddington ratio of the AGN outbursts fuelled by 
gravitational collapse. 



2.4 The Eddington ratio 

The Eddington ratios predicted by this model provide an- 
other method for comparison with numerical simulations 
and observations. Importantly, they also provide a way to 
check the self-consistency of the model, as outlined below. 

Broadly speaking, the form of energetic output from an 
AGN can be predicted from its Eddington ratio, defined by 



M Ed 



(10) 



where m is the black hole accretion rate. MEdd is the Ed- 
dington limited accretion rate determined by the balance 
between gravity and radiation pressure 

,-, -^Edd 47rGm p MBH 

J«Edd = 5- = , (J-J-) 

TjC r/CCTT 

where Z/Edd is the Eddington luminosity, Mbh is the black 
hole mass, or is the Thompson cross-section, m p is the pro- 
ton mass and the other symbols are as previously defined. 

By a nalogy with stellar m ass black holes in X-ray bina- 
ries (e.g. iKording et all [2008), accretion at rates less than 
~3% of the Eddington limit is radiatively inefficient so that 
the majority of the power output is in the form of kinetic- 
energy-dominated outflows of relativistic particles which are 
prominent radio synchrotron emitters. Outflows of this type 
are thought to couple strongly to the ambient gas. Con- 
versely, accretion above this critical rate is radiatively effi- 
cient, meaning that the power output is predominantly in 
the form of photons, rather than kinetic outflows of par- 
ticles. In this limit, radio jets may still be observed, but 
the efficiency of jet production is much lower than in th e 
radiatively inefficient regime (e.g. M accarone et alJ l2003h . 
Furthermore, in the radiatively efficient mode, it has been 
argued that only ~5% of the a ccretion power is availab le to 
heat t he su rrounding gas (e.g. ISiiacki fc Spr ingcl 200rj), see 
iKind (|2009T ) for a possible explanation. 

We note that the comparison with X-ray binaries is not 
exact, since it does not account for the existence of radio- 
loud quasars which produce powerful kinetic outflows at high 



© 2011 RAS, MNRAS 000, [JJ{8] 



AGN feedback triggered by gravitational collapse 5 



Eddington ratios. However, as explained below, we focus on 
low accretion rate objects so the distinction does not matter 
for the purposes of this model. 

To calculate the Eddington ratios, we use the well- 
known relation between the black hole mass and the stellar 
velocity dispersion: M BH » 1.3 x 10 8 (ct»/200 km s" 1 ) 4 M , 
where a* is the ce ntral stellar velocity dispersion (e.g. 
iTremaine et all 120021 ). Substituting for the black hole mass 
into equation (jTTJ yields the Eddington ratio for an AGN 
fuelled by gravitational collapse 



0.01 



io- 3 



200 kms- 



(12) 



where the assumed accretion efficiency is 77 ~ 0.1. Therefore, 
if /3/a < 10~ 3 , accretion following a gravitational collapse 
can be radiatively inefficient and produce kinetic outflows. 
Since a ~ 1, this implies that /3 < 10 -3 is required for kinetic 
outflow production. For larger values of /3/a the accretion 
will be radiatively efficient. 

The value of /3 is depends on the processes that govern 
how gas travels from kiloparsec scales onto the black hole. 
Since these processes are highly uncertain, we proceed by 
assigning /3 a single, empirical value that encapsulates the 
accretion physics in the full range of AGN environments. 
While this is a simplification, it ensures that the results are 
as transparent as possible. As indicated in equation (|12p . we 
find that f3 ~ 10 -3 permits a good agreement between the 



model and radio ob servations (|Best et al.l 120051 ; iBestl 120071 : 
IShabala et al.ll2008l ). Reassuringly, equation (|12p also shows 
that /3 ~ 10 -3 will lead to low Eddington ratios and, there- 
fore, kinetic outflows which produce radio emission. 



3 RESULTS 

From the assumptions and derivations outlined above, we 
present the expected duty cycles of radio AGN fuelled by 
gravitationally destablised gas, the corresponding duration 
of individual outbursts, and the associated heating rates re- 
sulting from shocks generated by the AGN outflow. 



3.1 AGN duty cycle 

In this section, we give a more general derivation for the 
AGN duty cycle. For a constant pressure cooling flow in 
which external heating is not important, the classical mass 
flow rate, gas temperature, and X-ray luminosity are related 
by 



7 ,y k B T 
fxrrip 



7-1 



(13) 



where fee is Boltzmann's constant and /jm p is the mean mass 
per particle. 

More generally, the net inflow of mass is determined by 
the difference between the time-averaged cooling and heat- 
ing rates. Observationally, this mass flow rate is estimated 
by fitting models to the X-ray spectra. Thus, we refer to 
the mass flow rate by its observational name, A/ spcc , but 
calculate it from the following 



H 



7 



■Ms, 



'/iffip 



(14) 



Rearranging equation (|14[) for the duty cycle, using equation 
0, gives 



. _ aGLx (. M spac 



(15) 



Consequently, in the limit that AGN heating exactly bal- 
ances gas cooling, M spcc /M c i as — ¥ 0, the duty cycle simplifies 
to 



7]f3c 2 a* (T 3 



(16) 



In terms of scaled quantities, the duty cycle can be expressed 
as 

-3 



8^0 



10 45 ergs" 1 J \200kms- 1 



(17) 



Equation (|17|l shows that AGN fuelled by the gravitational 
collapse of gas in more X-ray luminous clusters should ex- 
hibit larger duty cycles. That is, those AGN should be active 
for a gre ater proportion of time, as is inferred from obser- 
vations \Best et al.ll2007f ). However, the precise scaling of 
the AGN duty cycle in clusters is difficult to determine be- 
cause it is not clear exactly how the X-ray luminosity of the 
gas, and the central stellar velocity dispersion are related in 
clusters. Nevertheless, surveys indicate that the X-ray lumi- 
nosity of a cluster scales as Lx oc cr™ , where a c is the velocity 
dispersion of the cluster potential, and that n ~ 4 — 5 (e.g. 
iMahdavi fc Gelled |200iT ). Substituting the cluster L x - o c 
relation into equation (|17[) . gives 8 oc <7™/cr 3 . Therefore, if 
more massive BCGs are found in more massive clusters, the 
AGN duty cycle will tend to increase slowly with cr* , which 
would be in qu alitative agreement with observations (e.g. 
iBest et all 120071 ). Furthermore, assuming /3/a ~ 10" 3 , the 
normalisation of the duty cycle is en tirely consi s tent w ith 
the observational results presented bv lBest et al.l (|2007l ). 

In contrast, the scaling relations of field elliptical galax- 
ies yield a much simpler scaling of 8 with a». Under the 
assumption of an isothermal gravitational potential, the 
velocity dispersion is constant and independent of radius. 
Again, surveys indicate that the X-r ay luminosity scales as 
L x oc aT, where m ~ 8-10 (e.g. lO'Sullivan et all 120031 ; 
IMahdavi fe GellerilioOiT ). Using this fact, the AGN duty cy- 
cle would scale as 8 oc a™~ 3 ~ <r*. Since supermassive black 
hole mass scales as Mbh oc <t 4 , the duty cycle can then be 
said to increase as 8 oc Mrh , which i s consistent with ob - 
servational findings l|Best et al.ll2005l ; IShabala et al.ll2008l ). 
In addition, for a massive galaxy with an X-ray luminosity 
of 10 41 ergs _1 , and a* ~ 200kms _1 the duty cycle would 
be 8 « 1 0~ 4 (assum i ng /3/ a ~ 10~ 3 ), which i s also consis- 
tent with IBest et al l (|2005l ) and lShabala et al.l (|2008l ). Con- 
sequently, both the scaling and normalisation of the duty 
cycle are largely consistent with observations of radio AGN 
in field elliptical galaxies and BCGs. 



3.2 Number and duration of outbursts 

The AGN duty cycle is a useful quantity which provides 
information about how frequently an AGN is triggered in 
order for heating to balance the effects of radiative cooling. 
With additional assumptions, the duty cycle also provides 
some information about the duration of an individual heat- 
ing event, idyn, and its environmental dependence. If the 
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time between the onset of successive AGN outbursts is r, 
the number of heating events which occur during a galaxy 
lifetime, i agc , must he N — tagc/r = <5(t a ge/idyn)- Using the 
previous results, the expected number of AGN outbursts can 
be written 



N « 10 



10- 3 



i dyn / ^aj 

0.01 



10 45 



ergs" 



(18) 



200 kms 



By rearranging equation QB) , the duration of an outburst 



can be expressed in the form td y 



dtnge/N. Using equa- 



tion (|17l) for a cluster with Lx ~ 10 45 ergs _1 and <r* 
200kms _ , the typical outburst duration must be td yn ~ 
O.Oltagc Taking f ag e ~ 5 Gyr implies tdyn ~ 50Myr, which 
is sufficient to explain the fe atures AGN-blown bubbles ob- 
serve d in many clusters (e.g. lBirzan et a l. 2004; Dunn et al.l 
l2005h . 

The gravitational instability model of AGN fuelling 
also offers an explanation for the un expectedly large num- 
ber of compact ra dio sources (e.g. lO'Dea fc Bauml 1 19971 : 
IShabala et al.ll2008h that differs from the accre tion disk vari- 
ability model pro posed bv lCzernv et alj (|2009h . As shown by 
lAlexanderl (|2000h . the maximum stable length of an AGN 
jet propagating through an atmosphere depends on its power 
and the ambient density. For example, at constant jet power, 
a higher ambient density leads to a shorter stable jet length. 
Then, since the inflow of material due to gravitational insta- 
bility will significantly enhance the density in the vicinity of 
the black hole, this enhanced density may also plausibly con- 
fine and disrupt the resultant AGN jet on kiloparsec scales. 
If this is correct, the jet can only propagate further outwards 
once the gas has been sufficiently depleted by continued ac- 
cretion and star formation. 



4 DISCUSSION 

The model investigated above exhibits several features that 
are compatible with key observational characteristics of 
AGN feedback in both field elliptical galaxies and BCGs. 
Below we discuss some additional implications of the model 
which may be significant, but are harder to quantify. 



4.1 Star formation in the host galaxy 

Since only a very small fraction of the inflowing material 
reaches the black hole, the vast majority must either form 
stars, or be dragge d out of the galaxy by the AGN outflow 
itself. As shown in iPope et all i|201ol ) the mass of ambient 
material transported by an AGN-blown bubble is approxi- 
mately equal the mass of gas initially displaced by the bub- 
ble. In principle, this process can be extremely effective at 
removing material from the galaxy, dramatically reducing 
the quantity of gas available for forming stars. As a result, 
it is difficult to quantify star formation rates resulting from 
cooling-induced gravitational instability, because it is un- 
clear how much of the collapsed gas will be retained by the 
galaxy, and for how long. Our current best estimate is that 
the mass of material available for forming stars must fall be- 
tween two well-defined limits. The upper limit is the case in 



which all of the collapsed gas mass goes into forming stars; 
the lower limit is the case in which the vast majority of the 
collapsed gas mass is expelled by AGN feedback leaving no 
material available for forming stars. 

Given that the collapsing gas flows inwards on a dynam- 
ical timescale, tdyn, with a mass flow rate M, the collapsing 
mass can be written as M — Mtd yn - Then, using the fact 
that the instantaneous mass flow rate is M = al/(aG), the 
maximum mass of gas available for forming stars during each 
collapse, will be 



M = Mi 



dyn 



aG i 



10 1 



a * \ ( *dyn . , , 

200kms-V VK^yr/ Q ' 
(19) 

Since a fraction /3 is assumed to be accreted by the black 
hole, this leaves a maximum fraction of 1 — (3 available for 
forming stars, in the unlikely event that no material is re- 
moved due to feedback. 

However, we can estimate the mass of material removed 
by the AGN using simple considerations. The total energy, 
E, injected by an AGN as a result of accreting a mass /3M, 
will be E = rj/3Mc 2 . If this energy does work against the 
pressure, P, of the ambient gas, it will inflate a bubble with 
volume, V, given approximately by E = 4PV, assuming 
relativistic bubble contents. Using the definition of the gas 
pressure P — pksT/ (im p , where p is the ambient gas density 
and T is temperature, we can write E — iMdisksT/ fj,m p , 
where Mdi s = pV is the mass displaced in in flating the bub- 
ble. As previously noted, IPope et alj |2010) showed that a 
buoyant AGN-blown bubble will lift a mass comparable to 
Mdis/2 out from the central galaxy 0. 

The mass of gas lifted out from the centre of a galaxy 
by AGN feedback can be related to the mass of material 
that collapsed into the galaxy centre due to cooling-induced 
gravitational instability 



Mdi s /2 pm p c 2 



M 



Sk B T 



8 



P 



10- 



T 



10 7 K 



(20) 



According to equation (|20[) . an individual AGN outburst 
may remove more gas than actually flowed into the galaxy 
centre. While this may seem unphysical, it is not - feedback 
could remove all of the material that collapsed into galaxy 
and additional ambient gas. The fate of the gas lifted out 
of the galaxy will depend on the AGN power and the depth 
of the external gravitational potential. There are two main 
possibilities: 1) if the outflow injects sufficient energy, mate- 
rial will be permanently expelled from the galaxy; 2) if not, 
the outflow will temporarily lift material out of the galaxy, 
later allowing it to fall back inwards. 

2 It is difficult to estimate how high this material will be lifted 
because it depends strongly on the details of the bubble, the extra 
mass of mat erial it is carrying, and the properties of the ambient 
atmosphere. IPope et al.l ||2010|) demonstrated that a bubble lift- 
ing additional material will rise to a height at which the average 
density of the bubble, plus the lifted mass, is equal to the am- 
bient density. At this location the buoyancy force goes to zero, 
and the bubble cannot rise further unless it sheds some of the 
material. The greater the mass carried by the bubble, the less the 
bubble will rise. This effect directly limits the amount of energy 
extracted from the bubble. Because of this, it is not possible to 
say that the bubble will change the gravitational potential of the 
lifted mass by an amount defined by: E = (Mdi s /2) A<J>. 
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As shown in iPopd l|2009l ), AGN fuelled at a small frac- 
tion of the ambient cooling flow rate can power outflows 
that are capable of ejecting material from the potential of 
an elliptical galaxy. However, for the deeper gravitational 
potentials of galaxy groups and clusters (with virial tem- 
peratures greater than 1-2 keV) there are no black holes 
that are massive enough to sustain outflows that can expel 
gas from the potential. Despite this, typical AGN outflows in 
galaxy groups and clusters do still affect their environment 
by gently redistributing the gas within the gravitational po- 
tential. Furthermore, any material that is lifted out of the 
centre of the host galaxy will eventually fall back inwards, 
thereby becoming part of a fountain-like flow. 

Inspection of equation (|20[) leads to the expectation that 
AGN feedback should remove more ambient gas, per unit ac- 
creted mass, from cooler (lower mass) systems. As described 
above, this material can be permanently expelled from low 
mass elliptical galaxies. Thus, we conclude that star forma- 
tion in low mass elliptical galaxies would necessarily have to 
occur in a limited window of opportunity: after the gravi- 
tationally unstable material has collapsed into the centre of 
the galaxy, and before it has been removed by AGN feed- 
back. In other words, any reservoir of material available for 
forming stars will be short-lived, meaning that star forma- 
tion episodes are comparatively rare in such systems. 

When applied to galaxy clusters, equation (|20[) indi- 
cates that the mass of material lifted out of the host galaxy 
by AGN feedback will approximately equal the mass of ma- 
terial which collapses into the galaxy due to gravitational 
instability. As a result, AGN feedback may also be able to 
shut-off star formation in BCGs. However, the fountain-like 
flow described above would provide a reservoir of material in 
and around the galaxy which is available for forming stars. 
Consequently, is reasonable to expect some on-going star 
formation in BCGs. Interestingly, th is conclusion appears 
to be compatible with observations; iRaffertv" et al.1 (2006) 
found BCG star formation rates up to ~ 100 Mq yr _1 , which 
corresponds to 10 10 Mq over 10 s yr. 

Finally, considering the mass flow rates, we can show 
that the ratio of black hole and stellar mass growth rates 
must be M B h/M, = 13/(1 - P) « /9. This suggests that, 
as the age of the system becomes very large, the ratio of 
black hole mass to bulge mass should also tend towards (3. 
It is encouraging to no te that a value of (3 ~ 10 -3 agrees well 
with observations fe.g. lHaring fc Rbj2004l ') and is consistent 
with our earlier constraints obtained by compari ng the the- 
oretical AGN duty cycle w ith radio observations (|Best et al.l 
120051 ; IShabala et aTll2008l ) . 



4.2 



Possible implications for cosmological 
simulations 



The cooling-induced gravitational instability scenario for fu- 
elling AGN can only be captured by fluid simulations if the 
following are true: 1) the self-gravity of the gas is included in 
the calculations, and 2) the spatial resolution is fine enough 
to track the evolution of the gravitationally unstable region. 
For a collapsing mass M, with a velocity dispersion a,, the 
region of importance has a size 



Thus, for all plausible values of M and er*, the collapsing 
region is likely to be smaller than the resolved spatial scales 
in cosmological fluid simulations. We note that the density 
of these small, self-gravitating regions is likely to be high 



10" 22 gcm _3 for M 



10 10 M Q and a. ~ 200 km s" 1 ) 



R 



GM 



M 



10 10 Mq 



200 kms- 1 



kpc. (21) 



and scales as p ~ M /R 1, oc cr» /M 



5 SUMMARY 

The aim of this article has been to explore the properties of 
AGN feedback that is driven by the gravitational instability 
of hot gas in the locality of a supermassive black hole. While 
it remains unclear whether gravitational instability itself is 
responsible for delivering fuel to AGN, we have shown that a 
mechanism which behaves similarly could produce outcomes 
that are compatible with several key observations of radio 
AGN in massive galaxies and clusters. The main findings are 
summarised below: 

(i) Gas in the centre of a galaxy which periodically be- 
comes gravitationally unstable provides a way of linking the 
AGN fuelling rate to galaxy's large scale environment with- 
out requiring gas to flow 10s of kiloparsecs before reaching 
the black hole. 

(ii) According to this model, the AGN duty cycle scales 
as 5 oc Ly^/a'i, where Lx is the X-ray luminosity of the hot 
gas in a cluster and er* is the stellar velocity dispersion at 
the centre of the galaxy which hosts the AGN. 

(iii) Applying simple scaling relations to this model, we 
find that the duty cycle of radio AGN in massive galaxies 
should scale as 8 oc at oc Mbh, in reasonable agreement 
with observations. 

(iv) The region which collapses may be very small (< 
1 kpc) and is, therefore, difficult to capture directly in cos- 
mological numerical simulations. 

(v) The model predicts that the typical AGN outburst 
duration should scale as fd yn = St^/N oc Lx/(o"*-/V), where 
t age is the age of the cluster and N is the number of out- 
bursts during this time. Plausible values for t age and N imply 
that radio AGN outbursts fuelled by gravitational collapse 
can extend up to ~ 10 7 — 10 8 yrs, which is consistent with 
observations of AGN-blown bubbles. 

(vi) Using simple arguments we have shown that AGN 
feedback is likely to remove more gas, per unit accreted 
mass, from the centre of a lower mass galaxy than a BCG. 
By completely expelling material from lower mass elliptical 
galaxies, AGN feedback can dramatically reduce their star 
formation rates. In contrast, AGN feedback cannot com- 
pletely expel gas from a galaxy cluster. As a result, there 
may be more material in and around BCGs which is avail- 
able for forming stars. Thus, it may be more difficult to 
completely shut-off star formation in BCGs than in lower 
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